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ABSTRACT 

We present high speed photometry of the Sloan Digital Sky Survey cataclysmic vari- 
able SDSS J013701. 06-091234. 9 in quiescence and during its 2003 December superout- 
burst. The orbital modulation at 79.71 ± 0.01 min is double humped; the superhump 
period is 81.702 ±0.007 min. Towards the end of the outburst late superhumps with a 
period of 81.29 ± 0.01 min were observed. We argue that this is a system of very low 
mass transfer rate, and that it probably has a long outburst interval. 
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1 INTRODUCTION 

Cataclysmic variable stars (CVs) are identified in the Sloan 
Digital Sky Survey by spectroscopic classification of candi- 
dates selected on the basis of their colours, see Szkody et al. 
(2002). The discovery of SDSS J013701.06-091234.9 (here- 
after abbreviated SDSS0137) is announced in the second CV 
release of the Sloan Digital Sky Survey (Szkody et al. 2003) . 
The optical spectrum shows a blue continuum with broad 
absorption around narrow emission cores of H/3 and higher 
members of the Balmer series. There is also some indica- 
tion of broad TiO absorption features characteristic of an 
M dwarf primary (Szkody et al. 2003). These features in- 
dicate that it is a very low mass transfer rate (M) system 
where the accretion disc is so faint that most of the system's 
luminosity comes from the stars themselves. 

Szkody et al. (2003) obtained two hours of time re- 
solved spectroscopy and reported an orbital period (P or b) of 
84 min; the short P or (, is also consistent with a low M. The 
presence of TiO absorption, together with our photometry, 
indicates that this spectrum was taken during the normal 
quiescent state, rather than, say, during a decline from out- 
burst, where the hydrogen absorption may originate from 
the disc. 

We observed SDSS0137 at a quiescent V magnitude 
near 18.6 on seven nights, the system then brightened by 
~6 mag and strong superhumps were seen. This qualifies 



* E-mail: retha@mensa.ast.uct.ac.za 

f E-mail: pwoudt@circinus.ast.uct.ac.za 

| E-mail: warner@physci.uct.ac.za 

§ E-mail: gbolt@iinet.net.au 

% E-mail: jop@astro.columbia.edu 

II E-mail: armstrone@astro.columbia.edu 



SDSS0137 as an SU UMa type dwarf nova, see e.g. Warner 
(1995a) for a review of SU UMa stars. 

In Section [5] we present our observations in quiescence 
and outburst. These results are discussed in Section |3] and 
summarized in Section 0] 



2 OBSERVATIONS 

We obtained photometry of SDSS0137 at the Sutherland 
site of the South African Astronomical Observatory (SAAO) 
using the 1.0- and 1.9-m telescopes together with the Uni- 
versity of Cape Town charge coupled device (UCT CCD), 
O'Donoghue (1995) gives a description of the UCT CCD 
photometer and its reduction software. In addition, the 
SAAO CCD was used on three nights; unlike the UCT CCD, 
this is not a frame transfer CCD and hence there are a few 
seconds of deadtime between integrations. 

All the quiescent observations were made in white light; 
with the UCT CCD this gives photometry with an effective 
wavelength similar to Johnson V, but with a very broad 
bandpass. This, together with the nonstandard flux distri- 
bution of a CV, means that it is not possible to transform the 
observations onto any standard photometric system. By ob- 
serving hot subdwarf and white dwarf standards, we are able 
to give V magnitudes that are accurate to only ~0.1 mag. 

Most of the outburst observations were made with tele- 
scopes of the Center for Backyard Astrophysics (CBA), see 
Harvey et al. (1995). 



2.1 Quiescent light curves 

SDSS0137 was observed in quiescence on seven nights spread 
over almost four months, Tabled is a log of these observa- 
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Figure 1. The first four of the light curves of SDSS0137 in qui- 
escence, phased using equation^ Runs S7063, S7073, and S7088 
are shifted along the horizontal axis by 17, 53, and 169 cycles 
and down by 0.35, 0.60, and 1.00 mag respectively for display 
purposes. 



tions. We find a double humped orbital modulation, but no 
coherent signals on shorter time scales. The light curves are 
shown in Fig. Q Note that the relative strength of the fun- 
damental and first harmonic changes from night to night, in 
S7073, e.g., the first harmonic is only barely present. 

The lower scatter on the first three nights is caused 
in part by the fact that the SAAO CCD (used on these 
nights) gives a much larger field of view, and hence includes 
a larger number of comparison stars, than the UCT CCD, 
resulting in higher quality differential photometry. How- 
ever, at least some of the scatter in run S7088 is caused 
by intrinsic flickering. The UCT CCD's quantum efficiency 
peaks at A ~ 600 nm, but is still 25% at A ~ 350 nm 
(O'Donoghue 1995), whereas the SAAO CCD has very poor 
blue and ultraviolet (but high red and near infrared) sen- 
sitivity; see figure 17 of Woudt & Warner (2002). There- 
fore the UCT CCD observations are more receptive of short 
wavelength flickering. 

Fig. [5] is a Fourier transform of the first four runs, the 
fundamental and first harmonic give an unambiguous P or t = 
79.71 ± 0.01 min. The ephemeris for the times of maximum 



HJD max = 2452879.594 + d 055351(±7)£. 



(1) 



Fig. [3] shows the average light curve obtained by fold- 
ing runs S7060, S7063, S7073 and S7088 on this ephemeris 
after removing a linear trend from each of the light curves 
individually. 



2.2 Outburst observations 

In a system with M as low as indicated by the the spec- 
troscopic appearance of SDSS0137, one would expect to see 
only fairly infrequent dwarf nova outbursts. Nevertheless, af- 
ter leaving SDSS0137 at V = 18.6 mag on Dec. 17, we were 
lucky enough to find it at V = 12.5 mag four nights later. 
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Figure 2. The Fourier transform of runs S7060, S7063, S7073 
and S7088 combined. There is only one choice of alias that gives 
the correct fundamental and first harmonic relation, these peaks 
are marked by vertical bars at 0.2091 mHz (79.71 min) and 
0.4183 mHz (39.85 min). 




Figure 3. The mean quiescent light curve of SDSS0137 (plotted 
twice) using the first four runs folded on the ephemeris given in 
equation 1. 



Tablets a log of the observations during the outburst. This 
is the first recorded outburst. 

Smooth superhumps (with a peak to peak range of 
0.27 mag) were present on Dec. 21 already, implying that 
this was a superoutburst, that therefore SDSS0137 is an SU 
UMa type dwarf nova, and that the start of outburst was 
probably two or three days earlier. 

The outburst light curves are shown in Fig. [I] (a few 
noisy runs are omitted). Linear trends were removed from 
most of these light curves. Combining all the data from the 
first eight nights of the outburst gives the Fourier trans- 
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Figure 5. The Fourier transform of the first eight nights of the 
outburst. Vertical bars at 40.854 min and 27.231 min mark the 
first and second harmonic. 
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Table 1. Log of the observations in quiescence. 



Run No. 


Date (at the start 
of the night) 


HJD of first obser- 
vation (2452000 +) 


Length (h) 


tint (s) 


Telescope 


V (mag) 


S7060 


2003 Aug. 27 


879.54482 


3.20 


40 


1.0-m 


18.4* 


S7063 


2003 Aug. 28 


880.47545 


4.74 


40 


1.0-m 


18.4* 


S7073 


2003 Aug. 30 


882.51423 


2.51 


40 


1.0-m 


18.5* 


S7088 


2003 Sept. 5 


888.55917 


2.57 


40 


1.0-m 


18.5 


S7144 


2003 Oct. 5 


918.44016 


1.74 


.30 


1.9-m 


18.6 


S7168 


2003 Dec. 16 


990.28114 


2.32 


60 


1.0-m 


18.5 


S7172 


2003 Dec. 17 


991.27868 


2.27 


60 


1.0-m 


18.6: 



Notes: tint ls the integration time; ':' denotes an uncertain value; 'taken with the SAAO CCD. 
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Figure 4. The outburst light curves after linear trends were removed from some of them; most of these light curves were binned to 
improve signal to noise. The observations are phased according to the superhump period of 81.702 min. Individual runs are shifted along 
the horizontal axis by the number of cycles indicated, and arbitrarily displaced vertically. A few of the lower quality runs are not shown. 



form displayed in Fig. |H| The superhump period (Pah) is 
81.702±0.007 min; the first and second harmonic are also 
clearly detected. The light curves in Fig.0]are phased using 
this period. 

By the ninth night the superhumps were becoming no- 
ticeably phase shifted, another four nights later they had 
shifted by ~0.5 cycle and could be classified as late super- 
humps (e.g. Vogt 1983). The observed minus calculated (O- 
C) diagram (Fig. |SJ shows the transition from superhumps 
to late superhumps; there is no cycle count ambiguity in this 
diagram. 

Fig. shows the observations on HJD = 2453008 and 
one of the HJD = 2453002 light curves plotted on the same 
intensity scale. The amplitude and position of the late su- 
perhump maximum indicates that it did not develop out of 
anything visible in the superhump profile. Rather, the struc- 
ture in the superhump is its own repetitive profile. 

Fig. [SJis a Fourier transform of all the data taken after 
HJD = 2453004; it gives a late superhump period of 81.29 ± 
0.01 min. 



All the outburst light curves were searched for rapid 
quasi-coherent oscillations, but with no success. 



3 DISCUSSION 

It has often been noted that the standard model of CV evo- 
lution predicts the existence of a large number of short pe- 
riod, low M CVs, but that only very few such systems are 
known (e.g. Kolb 2002). Low M CVs are intrinsically faint, 
but many are being discovered in surveys that reach faint 
limits, such as the Sloan Digital Sky Survey. 

The orbital modulation of SDSS0137 is very similar to 
the double hump in WZ Sge (e.g. Patterson 1980). Three 
more examples of CVs with double humped orbital mod- 
ulations are WX Cet (Rogoziecki & Schwarzenberg-Czerny 
2001), SDSS J161033.64-010223.3 (Woudt & Warner 2004), 
and HS 2331+3905 (Araujo-Betancor et al. 2004). WZ Sge 
and WX Cet are both low M dwarf novae of long out- 
burst interval; no outburst has yet been observed in SDSS 
J161033.64-010223.3. 
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Table 2. Log of the outburst observations. 
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HJD 2452900.0 + 

Figure 6. O-C diagram computed by comparing the observed 
times of maximum with the times calculated from a sine function 
with a period of 81.702 min. This is consistent with a constant 
period up until HJD 2453003, after that the phase starts chang- 
ing. 
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Figure 7. Superhump and late superhump profiles plotted in 
intensity units. There is no vertical displacement between the 
two light curves, they are phased using the superhump period 
and the time of the first observed superhump maximum. The 
bottom light curve is transposed along the horizontal axis by 95 
superhump cycles. 
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Figure 8. The Fourier transform of all observations of the late 
superhump. Peaks corresponding to the fundamental and first 
harmonic arc marked. 



HJD of first 
observation 
(2452000 +) 


Length (h) 


k„t (s) 


Telescope 


V (mag) 


QQ 1 ^ 9Q039 


2.10 


5 


1.00-m 




QQ^ ^7788 
yyo.u / i oo 


4.93 


2 


1.30-m 


12.4 


nop. 09R74 




ID 


25-m 




QQfi 971 43 


0.96 


Q 


1.00-m 




QQfi R3391 


2.70 


5 


1.30-m 


12.6 


QQ7 fll 334 
yy / .uioo'i 


3.97 


15 


0.25-m 




QQ7 9Q90Q 
yy i .zy^uy 


3.76 


28 


0.36-m 




QQ7 4fi431 
yy i .'-tu i ±o J- 


3.62 


80 


0.25-m 


12.7 


QQ8 00849 


4.16 


15 


0.25-m 


12.8 


QQQ 01 974 
yyy.ui^ / *± 


3.81 


30 


0.25-m 


12.8 


1 000 0801 1 


2.21 


30 


0.25-m 




1 000 fil 81 9 


1.75 


45 


0.36-m 


13.0: 


1001.45554 


3.90 


80 


0.25-m 




1001.56973 


2.33 


60 


0.36-m 


13.0: 


1002.29981 


3.21 


28 


0.36-m 




1002.45077 


3.57 


80 


0.25-m 




1002.56369 


4.24 


5 


1.30-m 


13.1 


1003.27273 


1.93 


6 


1.00-m 


13.1 


1004.25491 


2.60 


28 


0.36-m 


13.2 


1006.23953 


0.88 


60 


0.41-m 




1007.08293 


1.63 


30 


0.25-m 


13.5 


1008.01126 


3.03 


30 


0.25-m 


13.6 


1011.01927 


3.13 


30 


0.25-m 




1011.30694 


0.87 


28 


0.36-m 


13.9 


1012.01266 


2.87 


60 


0.25-m 


14.7 



Notes: tint is the integration time; the magnitudes are nightly 
averages; ':' denotes an uncertain value; most observations 
were in white light. 



The double humped modulation in WZ Sge (and the 
other systems like it) has not been convincingly modeled, 
but is probably caused by the bright spot. An early expla- 
nation was that the bright spot has a large scale height and 
is thus visible above the inner accretion disc when it is on 
the far side of the disc (Krzeminski & Smak 1971; Robinson, 
Nather & Patterson 1978). Since then, it has been demon- 
strated that the disc in WZ Sge is optically thin (Skidmore 
et al. 2000; Mason et al. 2000) so that the bright spot can 
be seen through it. 

The second hump in SDSS0137 has lower amplitude 
during run S7073 than in any of our other runs. This has 
also been observed in WZ Sge (Krzeminski & Smak 1971); 
in WX Cet the second hump is sometimes not detected at 
all (Rogoziccki & Schwarzenberg-Czerny 2001). A possible 
explanation for this is that the optical thickness of the disc 
can temporarily increase so that the spot is (at least par- 
tially) obscured when it is on the side of the disc facing away 
from the observer (Krzeminski & Smak 1971). However, in 
the case of our run S7073, an increase in optical depth of the 
disc can't have been caused by an increase in M through the 
disc, since the system was 0.10 mag fainter than in run S7063 
and 0.14 mag fainter than in run S7060. Although there are a 
few features in the spectrum which may be attributed to the 
secondary, it does not contribute a large fraction of the total 
flux, therefore it seems very unlikely that the orbital modu- 
lation contains a substantial amount of ellipsoidal variation 
from the secondary. 
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Normal superhumps are thought to be caused by tidal 
stresses acting on an eccentric disc (e.g. Whitehurst 1988). 
Repetitive structure in superhump profiles similar to what 
can be seen in Fig. [7] and in more of the light curves in 
Fig. [I] has been observed in VW Hyi as well (Schoembs & 
Vogt 1980; Warner 1985). The superhump's fractional period 
excess is e = {P sh - P rb)/Porb = 0.0250 ± 0.0001. 

The interval between supermaximum and the first ap- 
pearance of superhumps is less than 4 days. This indicates 
that the outburst interval in SDSS0137 is probably years at 
most, rather than decades (Warner 1995b). 

A variation in bright spot brightness, resulting from the 
varying depth in the white dwarf potential well at which 
the stream impacts the disc, was initially proposed as a 
model for superhumps (Vogt 1982), but has since been used 
to explain late superhumps (e.g. Whitehurst 1988; Murray 
1996; Rolfe, Haswell & Patterson 2001), but see also Hess- 
man et al. (1992). Hellier (2001) proposes an alteration to 
the tidal thermal instability (Osaki 1989), namely that the 
tidal and thermal instabilities are uncoupled. The tidal ec- 
centricity can persist after the disc has come down from 
the thermal high state at the end of the outburst. The vis- 
cous dissipation - caused by tidal stresses - which results 
in superhumps, then decreases because the disc viscosity is 
now much lower. At this point the late superhumps should 
become detectable above the superoutburst light. It is there- 
fore perhaps surprising that the late superhump can not be 
seen in the HJD = 2453002 light curves. 

There are two other systems for which the transition 
from superhumps to late superhumps has been carefully 
tracked (VI 159 Ori (Patterson et al. 1995) and IY UMa 
(Patterson et al. 2000)). In both these examples O-C in- 
creases during the transition, in contrast to what is shown 
in Fig. E| 



4 SUMMARY 

We classify the CV SDSS0137 as an SU UMa type dwarf 
nova and report an orbital period of 79.71 min and a super- 
hump period of 81.702 min. SDSS0137 has a double humped 
orbital modulation, which is compatible with the very low 
M indicated by the spectral appearance. Late superhumps 
developed ~11 days after the start of the outburst. 
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